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We investigate the observational signatures of the holographic dark energy models 
in this paper, including both the original model and a model with an interaction 
term between the dark energy and dark matter. We first delineate the dynam- 

c3 . 

ical behavior of such models, especially whether they would have "Big Rip" for 
■ different parameters, then we use several recent observations, including 182 high- 

7— I \ 

quality type la supernovae data observed with the Hubble Space Telescope, the 

\Q • 

SNLS and ESSENCE surveys, 42 latest Chandra X-ray cluster gas mass fraction, 

y—i ' 27 high-redshift gamma-ray burst samples, the baryon acoustic oscillation measure- 

rs : 

ment from the Sloan Digital Sky Survey, and the CMB shift parameter from WMAP 



X 



three years result to give more reliable and tighter constraints on the holographic 



dark energy models. The results of our constraints for the holographic dark energy 
model without interaction is c = 0.748^QQQg, f2 m o = 0.276^Q'Q^g, and for model with 
interaction (c = 0.692l°;} 3 ^, O m0 = 0.28llg;^ ,a = -0.006+g;^, where a is an 
interacting parameter). As these models have more parameters than the ACDM 
model, we use the Bayesian evidence as a model selection criterion to make compar- 
ison. We found that the holographic dark energy models are mildly favored by the 
observations compared with the ACDM model. 
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I. INTRODUCTION 

It has been realized that our Universe is experiencing an accelerated expansion, as shown 
by several astronomical observations [l| • The acceleration of the Universe strongly indicates 
the existence of a mysterious exotic matter, namely the dark energy, which has large enough 
negative pressure and makes up the largest portion of the total matter in the current Uni- 
verse. The combined analysis of observational data suggests that the Universe is spatially 
flat, and consists of approximately 3/4 dark energy, 1/4 dust matter (cold dark matter plus 
baryons), and negligible amount of radiation 0,3,131. The simplest candidate of dark energy 
is the cosmological constant A (vacuum energy) which has the equation of state w = —1. 
The cosmological constant-cold dark matter model (ACDM) works very well, and is in agree- 
ment with a large number of resent observations. However, there are two problems in this 
scenario — the fine-tuning problem and the cosmic coincidence problem 4|. The fine-tuning 
problem asks why the vacuum energy density today is so small (10~ 47 GeV 4 ) compared with 
the theoretical value (10 74 GeV 4 ) from the quantum gravity. The cosmic coincidence prob- 
lem is that since the evolution of the energy densities of dark matter and dark energy are 
so different during the expansion of the Universe, why are they nearly equal to each other 
today? 

The dark energy may be a problem which has to be solved in the context of quantum 
gravity [5]. In the classical gravity theory, the dark energy density (cosmological constant) 
can be an arbitrary value. However, a complete theory of quantum gravity should be capable 
of determining the properties of dark energy such as the energy density and the equation 



of state |5[. The holographic dark energy model is an attempt to app 
principle of quantum gravity theory to the dark energy problem 
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the holographic 



It is well known that the holographic principle is an important result from the explo- 
rations of the quantum gravity theory and string theory [lOj], enlightened by investigations 
on the properties of black holes. For an effective field theory in a box of size L with UV 
cut-off A c , the entropy S scales extensively as S ~ L 3 A^. However, considering the pecu- 
liar thermodynamics of black hole [lj| , the maximum entropy in a box of volume L 3 may 
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behave nonextensively, growing only as the surface area of the box, i.e.S < Sbh = nM 2 x L 2 
(Bekenstein entropy bound). This nonextensive scaling suggests that quantum field theory 
breaks down in large volume. To reconcile this breakdown with the success of local quan- 
tum field theory, Cohen et al. 6j proposed a more restrictive energy bound. They pointed 
out that in quantum field theory a short distance (UV) cut-off is related to a long distance 
(IR) cut-off due to the limit set by forming a black hole. In other words, if the quantum 
zero-point energy density p\ is relevant to a UV cut-off, the total energy of the whole sys- 
tem with size L should not exceed the mass of a black hole of the same size, thus we have 
L 3 p\ < LM 2 V When we take the whole Universe into account, the vacuum energy related 
to this holographic principle lOj is viewed as dark energy, usually dubbed holographic dark 
energy. The largest IR cut-off L is chosen by saturating the inequality so that we get the 
holographic dark energy density 



p dc = 3c 2 M^L- 



(1) 



where c is a numerical constant, and M p \ = l/v8vrG is the reduced Planck mass. If we 
take L as the size of the current Universe, for instance the Hubble scale H~ l , then the dark 
energy density will be close to the observational result. However, Hsu jg] pointed out that 
this yields a wrong equation of state for dark energy. Li 9] subsequently proposed that the 
IR cut-off L should be taken as the size of the future event horizon 

dt' r°° da' 



RphiCl) 



(2) 



, a(f) "./„ Ha' 2 ' 

Then this problem can be solved, and the holographic model can thus be constructed sue 
cessfully. For extensive studies of this model (HDE), see Ref. 
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3, H, Q, 0, Q, 



23]. The holographic dark energy may also have interaction with matter, 
this is the so called interacting dark energy model (IHDE) H Q. Both the HDE and 



HDE model can be tested by cosmological observations, such as the type la supernovae 



24j , CMB 25|, |26|, |27| , X-ray gas mass fraction in the clusters [28j] , differential ages of pas - 



sively evolving galaxies [29J and combinations of SN la, CMB and LSS data 
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33|. 



32j, namely c = 0.85±° 



+0.18 
02 



Currently, the tightest constraint on the HDE model is given by 
and fl m0 = 0.271q;q3, but c > 1 is still allowed under this constraint. In addition, by us- 
ing the x 2 statistic as a model comparison technique, Ref.[32| suggests that HDE model is 
equally favored by the current observational data compared with the ACDM model. 
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In this paper we consider the dynamical behavior of the holographic dark energy models, 
particularly what is the condition for the Universe to end in the so called "Big Rip" \s4\. 
we then use several new data sets to constraint the holographic dark energy models. These 
includes a sample of 182 high quality SN la data (Despite coincidence in number, this is 
not exactly the same data set used by Gold06 [35(), a sample of 42 latest X-rays gas mass 
fraction (f gas ) data [38J|, a sample of 27 gamma-ray burst (GRB) data 39] generated from one 
of the tightest correlations-the - E eolation Q, the baryon acoustic oscillation 

measurement from the Sloan Digital Sky Survey [411], and the CMB shift parameter from 
WMAP 3 years result. With these more diverse set of observations, we can obtain constraints 
which are not only tighter, but also more reliable. 

We compare the observational fit of the HDE and IHDE models with the ACDM model. 
In previous works, the y 2 statistic has been used for comparison. However, as the holographic 
dark energy models (with or without interaction) have more parameters than the ACDM, 
the x 2 may not be a suitable criterion for making comparison. Here we use the Bayesian 
Evidence (BE) as the information criterion to assess the strength of the holographic models. 

The paper is organized as follows: In §2, we review our models of holographic dark 
energy, and derive the evolution equations. We also discuss qualitatively the behavior of the 
HDE and IHDE models, particularly the impact of the interacting term on the fate of the 
Universe. In §3, we present our method of analysis and the data used (§3A), then briefly 
discuss the Markov Chain Monte Carlo (MCMC) techniques (§3B) and the model selection 
criteria (§3C). In §4, we show the results of our constraints, and using the Bayesian evidence 
we make comparisons between the holographic dark energy models (including the HDE and 
IHDE) and the ACDM model. We summarize our result and conclude in §5. 



II. THE MODELS 



For a spatially flat (the flat geometry is assumed throughout this paper) Friedmann- 
Robert son- Walker (FRW) Universe filled with matter component p m and holographic dark 
energy p de , the Friedmann equation reads 

3M 2 # 2 = p m + Pde . (3) 

where pde = 3c 2 MpjL~ 2 , L is the future event horizon, i.e. L(t) = R eh (a). 
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When we consider the interaction team between the two dark components, the conserva- 
tion equations can be written as 

p m + 3Hp m = Q, (4) 

p de + 3if (1 + w dc )p dc = -Q, (5) 

where w de is the equation of state of holographic dark energy. The form of the interaction 
Q is not unique. Here, we consider an interaction term of the following form (similar to but 
slightly different from those given in jl?!]): 



Q = 3aHp dci 



(6) 



This form of interaction has been discussed in some literatures on inflation and reheating. 
For example, in the warm inflationary model 42|, in which the scalar field's energy is 
transferred to the matter due to scalar field oscillations, there is an interaction term which 
is throughout the inflationary regime (not just after slow-roll), so that the energy of the 
scalar field is transferred to the matter content continuously and the matter content is not 



driven to zero 43]. Also, in string theory, a similar interaction term arises in the Einstein 
frame which depends on the dark energy density. Moreover, a term of the form ([6]) could 
be motivated by analogy with dissipation, for instance, a fluid with bulk viscosity may give 

:hat 



rise to a term of this form in the conservation equation [43l . |44|. We note, however, t 



interactions of other form is also possible and has been discussed in the literature 45 



46 



47| 



Using the definition of holographic dark energy (CQ) and the relationship pa e = 3H 2 M^Qde 
[in which f2d e is the fractional dark energy density), taking the derivative with respect to 
x = In a, we obtain 

dpde c i\ /r2 rr2o 1 '1 V^de\ 



Pdc 



-QM^H 2 n dc (l 



de- 



dx *" C 

Considering the derivative relationship between t and z: 4 = H-4- = —H(l + z)4- and Eq. 
(JHJ) and (jHJ), we have the following equation 



Pde + 3{1 + w de )p dc = -3ap de . 



(8) 



After taking derivative about p de = 3H 2 M^Q dc and substitute Eq.([7j) into it, we obtain 

n: 



H 

IT 



2tt 



,|,. V^dT 



1. 



do 



C 



(9) 
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On the other hand, using the Friedmann equation H = — -^hr^P + p)(p and p are the total 

pi 

energy density and pressure), and substitute H = H H and w de from Eq. (E]) into it, we 
could get 

H' 1 3 1 1 3 

= -^c - - + -fife + -"fide- (10) 

Combining Eq. (Q and ([TO]) , we find the differential equation for f^e 

d ^^L = n dc [(i - n dc )(i + -y/HZ) - 3afi de ], (ii) 
dx c 

i.e. 

+ ^[(1 - fi dc )(l + -V^Z) - 3afi dc ] = 0. (12) 
dz 1 + z c 

Consequently, the differential equation for Hubble parameter H(z) could be written as 

dH H(z) A 



dz 1 + 2 [^(l + 3 tt + -v^)--]. (13) 

Eq. (fl2l) and (FTBl) can be solved numerically to obtain the expansion rate H(z). 

We now study the behavior of the holographic dark energy models. The evolution of dark 
energy can be understood by inspecting its effective equation of state, for which we have 

p' de + 3{l + w eS )p de = 0. (14) 

Using Eqs.([7D and (JSJ) we obtain 



f \ i 12 y/Qdc 

w eS {z) = w de + a = -- - . (15) 

o o C 

It is interesting to note here that a (i.e. the interaction rate of dark energy and dark matter) 
does NOT show up explicitly in this equation, but only implicitly affects the result through 
its effect on Q dc , so this equation applies to both HDE and IHDE models. For the flat 
Universe with dark energy, the dark energy would eventually dominate the density, Q dc — > 1, 
hence the effective equation of state of the holographic dark energy evolves dynamically, and 

U> e ff — > — - — — . (16) 
3 3c v ; 

It is obvious from Eqs. (ITB]) - (IT6T) that for the HDE model, the equation of state depends 
on the parameter c. For c > 1, the equation of state w > — 1, and the holographic dark 
energy behaves as quintessence. For c < 1, w < — 1 could be realized. For such models, the 
Universe would end in a Big Rip. 
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FIG. 1: Equation of state for the selected values of (c, a) or the illustration purpose. Here we set 
fimO = 0.27. 

For IHDE models, we illustrate the behavior of the model by choosing some representative 
values of a and c, and plot the evolution of equation of state w(z) in Fig. [TJ On the left of 
Fig. [JJ we show the evolution of equation of state w as a function of z for fixed value of c 
(c = 1.05) and different values of a. Initially, all models have —1 < w e ff < 0. If a = 0, the 
model is reduced to the case of HDE, for which the equation of state w — ■> — 1 (cosmo logical 
constant) as z — > — 1 (a — > oo). For a > 0, energy is transferred from dark energy to dark 
matter, making the effective equation of state of the dark energy greater than —1, so the 
dark energy behaves as quintessence-like. For a < 0, energy is transferred from dark matter 
to dark energy, so the effective equation of state of dark energy crosses —1, i.e. it exhibits 
quintom behaviors. On the right side of Fig. [TJ we plot the equation of state for fixed a 
with different values of c (a = —0.11 in this particular case). As c increases, the equation of 
state depart from —1 and increases. For large value of c the quintom divide is not crossed, 
and the model again behaves like a quintessence. In Fig. [2] we plot pdc{z)- On the left of 
Fig. El we plot the special case of HDE (a = 0). we find that for < c < 1, the equation 
of state will cross —1 and behave like a "quintom", so that phde diverges in finite time; for 
c > 1, the equation of state is greater than —1, and the dark energy density decrease as 
time passes. On the right of Fig. [21 we plot behavior of the system with different values of 
a with c = 1. For sufficiently negative a, the IHDE density diverges as energy is transferred 
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FIG. 2: Energy density of IHDE and HDE model. Left: Energy density of HDE . Right: Energy 
density of IHDE. Here we set fi m o = 0.27. 

from the dark matter to dark energy, and eventually we would run into the "Big Rip" . 



III. METHODOLOGY 
A. Data analysis 

We utilize several data sets to constrain the parameters of the holographic dark energy 
model, including a selection of 182 high-quality type la supernovae, the baryon acoustic 
oscillation measurement from the Sloan Digital Sky Survey, 42 latest X-ray gas mass fraction 
data from Chandra observations, 27 GRB samples generated with E pea k — E 1 correlation, 
and the CMB shift parameter from WMAP 3 years result. 



1. Selected SN la data set 



For t 
sample 



le SN la data, our sample includes 77 highest quality data from the Riess Gold 06 



351 ] (30 HST supernovae and 47 high quality SNLS supernovae 



36j)[48|, and 105 



ESSENCE supernovae [33] (60 observed by the ESSENCE team and 45 nearby supernovae 
re-analyzed by the ESSENCE team), the total number is 182. We use the light curve fitters 



MLCS2k2 



49 



501 ] (the other light curve fitter SALT [51] give results which are consistent 
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with MLCS2k2 35|, l37j), this algorithm avoids the need of normalization 52|. The redshift 
of this data set reaches 1.755. 

The likelihood function of the parameters can be determined from \ 2 statistics, and for 
SN la data 



182 



i=l 



(17) 



where the theoretical value of distance modulus fith{z) is given by 

Hth{z) = 5log l0 d L {z) + 25 

= 5 log 10 D L {z) - 5 log 10 ho + 42.38, 



and 



D L (z) = — x d L (z). 
c 



Here Hq = 100/iokms 1 Mpc 1 and the luminosity distance cLl can be written as 

cdz' 



d T ,(z) 



1 + z) 



H(z' 



(18) 



(19) 



(20) 



2. BA O measurement from SDSS 



The baryon acoustic oscillation signatures in the large-scale clustering of galaxies can 
be seen as a standard ruler providing another way to explore the expansion history of the 
Universe. We use the measurement of the BAO peak from a spectroscopic sample of 46,748 
luminous red galaxies (LRGs) observations of SDSS to test cosmology [4l|, which gives the 
value of A = 0.469(n s /0.98)-°- 35 ±0.017 at z BAO = 0.35 where n s = 0.95 [53]. The expression 
of A can be written as 

1 2/3 



A 



mO 



(H(z BAO )/H y/ s 



^BAO 



dz' 



zbao Jo H(z')/H 



(21) 



and the Xbao i s 



, A- 0.469K/0.98)- - 35 
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3. Latest X-ray gas mass fraction data from Chandra 



The X-ray gas mass fraction in the largest, X-ray luminous, dynamically relaxed clusters 
of galaxies provides a fair sample of the matter content of the Universe. It could give a 
constraint to the geometry of the Universe with the relation f gas cx dj 5 , under the assumption 



33 



that this fraction should be approximately constant with redshift 
the latest da ta derived from 42 relaxed clusters by Allen et. al 
the redshift of this sample ranges from 0.05 to 1.1. 
% the X/ gas is given by 



54 



551 ] . Here we use 
, from Chandra 



observations 



38 



Following Allen et. al [38 

[ xACDM 




(23) 



and the model fitted to the reference ACDM (Q m o = 0.3, Q 



AO 



ACDM/ \ = fv^T^Vf 

gas 1 ' l + s(z) 





~di CBM (z)~ 




d A (z) 



0.7) data is 

1.5 



(24) 



where dA,(z) and d\ CT>M (z) are the angular diameter distances to the clusters in the test 
model and reference model, 

cdz' 



w-VTTd-m- (25) 

The parameter b(z) = bo(l + a&z) in Eq. (124"1) is the bias factor which parameterizes 
the redshift-dependent deviation of the baryon fraction measured at r25oo from the Universe 
mean with 0.65 < bo < 1.0, —0.1 < a;& < 0.1; The factor s(z) = so(l + ot s z) models the 
baryonic mass fraction in stars, and sq = (0.16±0.05)/?-75 5 , —0.2 < a s < 0.2. The A accounts 
for the change in angle subtended by r 2 5oo as the reference cosmological model is varied: 



A=r- 



n ACDM 



(26) 



\ #2500 / V [H(z)d A (z)] J 

and n is the slope of the / gas in the region of r 2 5oo measured in the reference ACDM model, 
which takes the value r] = 0.214 ± 0.022. The parameter 7 represents the effect of non- 
thermal pressure support in the clusters, which ranges from 1.0 to 1.1; The factor K is a 
'calibration' constant which accounts for residual uncertainty in the accuracy of the instru- 
ment calibration and X-ray modeling, and we take K — 1 ± 0.1. 
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4- Gamma-ray bursts data 



57 



he GRB data may be a good complement to the other observational data 39j, |40j, |56, 



5§J, 



59 



60 



61] , such as SN la data. They have very large redshift distribution and can 
be observed at much higher redshift, thus provide an effective way to detected the evolution 
of the dark energy. 



Our GRB data set is constituted of 27 GRB samples in Ref. [39j. They are generated from 
the Ep ea k — Ky correlation (Ghirlanda relation [401]). which is one of the tightest correlations 
for GRB. The redshift of this data set reaches 6.29. 
The Xgrb takes the form of: 



Xgrb . , 



in which a distance modulus fi bs(z) estimated from the observational data can be calculated 

as 

fi obs (z) = 5 log 10 (d Lo J + 25, (28) 
with a estimated luminosity distance d Eoba expressed in unit of megaparsecs, and 

d Lobs = [£ 7 (1 + z)/(AnF beam S bolo )} 1/2 . (29) 

Here F beam is the beaming factor, S bo \ is the bolometric fluence of the burst, and the 
collimation corrected energy E 1 can be fitted by 

log£ 7 = a + b\og[E peak (l + z)/300A;eV], (30) 

where a = 50.57 and b = 1.63. Its uncertainty is 

*l g E, = ol + fa hg[E peak (l + z)/300keV}) 2 + (0A3A3ba Epea J E peak ) 2 + a 2 E ^ sysJ (31) 

where the la uncertainties in the intercept and slope are a a = 0.09 and o b = 0.03, and the 
best estimated cr_g S2/s is 0.16. The <Ji in Eq. ( |271) can be estimated as 

<n = [(2.5a logi?7 ) 2 + (l.086a Sbol JS bolo ) 2 + (1.086^ eom /F 6earo ) 2 ] 1/2 - (32) 
At last, the expression of Hth{z) m Eq. (1271) is given by Eq. (TI8"]) . 
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5. CMB shift parameter from WMAP 3 years result 

The shift parameter is a very good complement to the previous data set because the very 
large redshift distribution (zcmb = 1089) can reflect the evolution of the dark energy. The 
shift parameter R is derived from the CMB data takes the form as 

/* Z CMB A z ' 

R ^^W„ (33) 



The WMAP3 data gives R = 1.70 ± 0.03 [62} , thus we have 



2 /R-1.70\2 
Xcmb = • (34) 

To break the degeneracy and explore the power and differences of the constraints 
for these data sets, we use them in several combinations to perform our fitting: 
SN sel + BAO, SN scl + BAO + f gas , and SN sel + BAO + f gas + GRB + CMB. 

B. MCMC 

Markov Chain Monte Carlo (MCMC) techniques are widely used to generate random 
samples to simulate the posterior probability of the parameters given the data sets. This 
method has several advantages over grid-based approach. Most importantly, the computa- 
tional time cost increases approximately linearly with the number of parameters, so even for 
a lar ge n umber of p arameters the estimate can be done within an acceptable computation 



time 



63 



64j, 



65 



66). 



We use the Metropolis-Hastings algorithm with uniform prior to decide whether to accept 
a new point into the chain by an acceptance probability: 



a(0 n+ i|0 n ) = min 



p(0 n+ i|d) q(0 n |0 : 



n+lj 



p(0„|d) q(0„+i|0 E 



where p(9) is the prior probability distribution and q(# n +i|#n) is the proposal density of 
proposing a new point 8 n +i given a current point 9 n in the chain. If a = 1, the new point 
9 n+ i is accepted; otherwise, the new point is accepted with probability a. The trials are 
repeated until a new point is accepted, and then we set 9 n = 9 n+ \. In our computation, 
we set a Gaussian-distributed proposal density for every point which is independent of the 
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position on the chain, so that q(# n +i|#n) and q(9 n \9 n+1 ) are canceled, and consider the 
uniform prior and Bayes' theorem we get 

aMHmin J. ,35) 

Here C(d\9) is the likelihood to obtain the data set d given the parameter set 8, and usually 
can be written as 

£(d\0) = -T^e-^ 2 . (36) 
For our three data sets, the y 2 are 



x 2 



Xsn sc1 + Xbao 



Xsn s „ + Xbao + Xf gas ( 37 ) 



XsN sel + Xbao + xL s + Xgrb + Xcmbi 



respectively. We assume uniform prior for the parameters of our models within the given 
ranges as following: Q m0 G (0,1), c G (0,2), a G (-0.2,0.2) and h G (0.4,0.9). In 
particular, when we use f gas data, the parameters coming from this data set 6>d a ta = 
{fife, s , a si b , at,, r], 7, K }, are also included in our MCMC fitting process. 

The thermodynamical bound c > \JVt£> is NOT assumed a prior, as we wish to assess the 
value of the data fairly. 

We generate six chains for each case we study, and about one hundred thousand points 

n 

are sampled in each chain. The form of proposal density we use is described in Ref. [671 ] . 
After the convergence determined by Gelman and Rubin 68( criterion and thinning the 
chains, we merge them into one chain which consists of about 10,000 points used to simulate 
the probability distribution of the parameters. 



C. Model Comparison 

For comparing different models, one must choose a statistical variable, the x 2 min is the 
simplest one and is widely used. However, for models with different number of parameters, 
the comparison using x 2 may not be fair, as one would expect that models with more 
parameters tends to have lower x 2 - 



The Akaike information criterion (AIC) 69] 



AIC = -2 1n£ max + 2A ; (38) 
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includes the penalization of the number of parameters, where C max is the maximum like- 



lihood and k is the number of parameters 



70 



71] . However, the size of the data is not 
embodied in the AIC, if there is a large number of data, then the reduction in x 2 due to 
the additional parameters may be very large and the 2k term in the AIC smaller could not 
compensate it [72]. 



The Bayesian information criterion(BIC) 73] can be written as 



BIC = -2 lnC max + k \nN, 



(39) 



where N is the number of data. The B 



much if given large number of data 



0, 



C tends to penalize the number of parameters too 



72, 



74). 



We will use the Bayesian evidence (BE) as a model selection criterion. The Bayesian 
evidence of a model M takes the form 



BE = / C(d\9,M)p(9\M)d9, 



(40) 

where C(d\8, M) is the likelihood function given the model M and parameters 9, and p(8\M) 
is the priors of parameters. The BE may be the best model selection criterion, as it is the 
average of likelihood of a model over its prior of the parameter space and automatically 
includes the penalties of the number of parameters and data, so it is more direct, reasonable 

AIC and BIC in model selection [70 



and unambiguous than the x 



2 

mini 



75 



76|, 



77 



78 



79 



80, 



8l| (For a connection between BE and Xmin analysis, see Ref. |s2]). The logarithm of BE 
can be used as a guide for model comparison (Jeffreys 1961), and we choose the ACDM as 
the reference model: A In BE = \nBE mo d e i — IuBEacdm- The strength of the evidence for 
the model is considered according to the numerical value of BE: 



A In BE < 1 
1 < A In BE < 2.5 
2.5 < A In BE < 5 
A In BE > 5 



Weak 
Significant 

Strong to very strong 
Decisive 



(41) 



We use the nested sampling technique to compute BE [7_8j, l83 |. 
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IV. RESULTS 

A. The HDE Model 

In Table I, we give the best fit value and the la error of the HDE model parameters, as 
well as the value of In ABE for the three data set combinations. We plot the probability 
distribution function (PDF) of the HDE model in Fig. [3j We can see that the best fit of Q m o 
are almost the same (around 0.27) for all data sets. The best fit of c varies slightly across 
the different data sets, it is 0.761 for the SN+BAO data set, but decreases slightly when the 
/ gas , GRB and CMB data are included. However, for all data sets, we have c < 1 at more 
than 1.5cr (see Fig. [3]). The PDF of the parameter distribution is smoothly distributed. 




0.25 0.3 0.35 0.4 0.6 0.8 1 1.2 



FIG. 3: Probability distribution function (PDF) of the parameters for the fits of HDE model. Left: 
The PDF for parameter 0, m Q. Right: The PDF for parameter c. 



TABLE I: The fitting result for the HDE model 





SN + BAO 


SN + BAO + f gas 


SN + BAO + f gas + GRB + CMB 


C 

A In BE 


n 97 o+0.020 

n 7fii+ - 154 

u - ' Di -0.117 

0.09 ±0.12 


2 70+ - 021 
u - z ' u -0.018 

n 7^+0.130 

u -' 4O -0.101 

0.63 ±0.18 


0.276±g3JS 
0.65 ±0.18 
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FIG. 4: The contour maps of c vs. O m0 for HDE with 1<t(68.3%), 2cr(95.5%) and 3cr(99.7%) 
confidence levels. 

The parameter c plays an essential role in determining the evolution of the HDE model. If 
c = 1, the dark energy equation of states would asymptote to that of a cosmological constant 
and the Universe would enter the de Sitter phase in the future; if c > 1, the equation of 
state of dark energy would always be greater than —1, it would behave as quintessence dark 
energy; if c < 1, initially the equation of state of HDE would be greater than -1, but it 
would decrease and eventually cross the "phantom divide line" (w = —1) as the Universe 



expands, acting clS db quintom 



84] 



We plot the contour maps of cvs. f2 m0 for HDE with la(68.3%), 2a(95.5%) and 3cr(99.7%) 



confidence levels in Fig. HI Our constraint is tighter than previous ones, e.g. Ref. [30j], as 
we have used more precise data in our fitting. The center of the best fit is located at c < 1, 
but there is still a fair fraction of allowed parameter space in which c > 1. The evolution of 
w and pde in HDE models with the best fit parameters for the three data set combinations 
are shown in Fig. [5j For these cases, as we expected, the dark energy diverges in finite time 
and the Universe ends with a Big Rip. 

Moreover, for all the four data sets, we have c < 1.2 at more than 3a (see Fig. [3] and 
Fig. Hj), which is rather consistent with the possible theoretical limit of parameter c from 



the weak gravity conjecture (see |23j)- 
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The HDE model fits about equally well (In BE=0.09) as the ACDM when we only use 
the SNIa and BAO data. With f gas , GRB and CMB data added, it fits mildly better than 
the ACDM, but with the data presently available the difference is not significant (In BE 
=0.63 ~ 0.65). 
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FIG. 5: Evolution of the dark energy in the best fit HDE model. Left: Equation of state. Right: 
Dark energy density. 



B. The IHDE Model 



TABLE II: The Fitting results the IHDE model 





SN + BAO 


SN + BAO + f gas 


SN + BAO + f gas + GRB + CMB 


C 

a 


n 979 +0.023 
u - z ' z -0.022 

0.592j£™ 
-0.020i°;^I 


n 97 c:+0.021 
u - z '°-0.021 

U - DD ' -0.164 

0.0681°;?! 


o-28i±8:S? 

0-692±8;ig| 
-n nnfi+°- 021 

U.UUD_ 024 


A In BE 


0.41 ±0.12 


0.70 ±0.18 


0.75 ±0.18 



In Table II, we give the best fit value and the la error of the IHDE model parameters, 
as well as the value of In ABE for the three data set combinations. We plot the probability 
distribution function (PDF) of the IHDE model in Fig. El We also plot the Q m o — c contours 
in Fig. [3 The evolution of the effective equation of state and the relative density of the 
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FIG. 6: The PDF of the parameters for the fits of IHDE model. Left: The PDF for parameter c. 
Right: The PDF for parameter a. 

dark energy are plotted in Fig. [HI Similar to the case of HDE models, the best fit of Q m o 
for the IHDE models is around 0.27 0.28 for all three data set combinations, as can be seen 
from Fig. [7J However, for the different data set combinations, the distribution of c and a 
are fairly different (see Fig. [6] and Table II). The peaks of the PDF for a are different for 
the three data set combinations. Furthermore, the evolution of the equation of state for the 
three data set combinations are also very different (see Fig. [HJ. 

Does this mean that the three different data sets indicate very different expansion behav- 
ior? If so, this would indicate that the three data sets may be inconsistent with each other. 
However, we have seen that for the case of HDE the three different data set combinations 
yield similar fitting parameters, indicating that they are basically consistent with each other. 

To understand the origin of this difference, we plot the expansion rate in Fig. [9] for the 
best fit model parameters of the three data set combinations. While there are differences, 
one can see qualitatively the three set of curves are similar to each other. This shows that 
the different data sets are not inconsistent. 

The reason of the difference seems to be due to parameter degeneracy in the IHDE model. 
In this model, we have two parameters c and a, it appears that different combination of 
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0.5 1 1.5 Z 



FIG. 7: The contour maps of c vs. ft m0 for IHDE with 1<t(68.3%), 2ct(95.5%) and 3ct(99.7%) 
confidence levels. 

these two parameters may lead to similar dynamical behavior. At first sight, this seems to 
be in disagreement with our analysis in §2, where the effective equation of state w e // of the 
dark energy was given in terms of these two parameters. However, the equation of state 
w e ff of the dark energy is not the only thing affecting the expansion. In this model, the 
interaction between dark matter and dark energy induces a non-zero effective equation of 
state for the dark matter. Thus, although the effective equation of state for the dark energy 
looks very different, the change in the dark matter equation of state compensates part of this 





-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



z z 



FIG. 8: Effective equation of state and density evolution for the best fit IHDE models. 
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FIG. 9: The expansion rate H(z) for the best fit parameters of the IHDE model. 

difference. Further, this model still has some shortcomings since fide could quantitatively 
evolve to region larger than 1 if a < 0, i.e. if dark matter decays into dark energy, because 
the interaction term does not concern the energy density of dust matter [4^]. At this case, 




FIG. 10: The contour maps of a vs.c for IHDE with 1<t(68.3%), 2a(95.5%) and 3cr(99.7%) confi- 
dence levels. The black dot-dashed curve denotes w e R = — 1 when z — ► — 1 with fiiHDEO = 0.73, 
and the region below (over) it means w e g will (not) cross -1 during finite time. 
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the IHDE model is only effective in the period prior to the Big Rip. We will improve our 
work and discuss the possible interaction forms in the next several papers. 

The contour map for c and a is plotted in Fig. [TQ1 From the Fig, the contours correspond- 
ing to the data sets SN+BAO+/ gas +GRB+CMB are much tighter than the other data sets. 
This is because we use the data set GRB and CMB, which has very large redshift distribu- 
tion so they break the degeneracy of the parameter c and a. We also mark the a — c values 
for which w e ff = —1 as a dashed line. This forms the dividing line between quintessence-like 
and phantom-like behavior. For the best fit parameters of all three data set combinations, 
c ~ 0.6 < 1. The value of a varies more, but all consistent with being within 1.5<r, there is 
no strong evidence for the presence of interaction. For the SN+BAO+/ gas -|-GRB+CMB, the 
PDF of parameter a and the best fit values strongly suggest the evidence for the interaction 
is very weak. In any case, for all three data set combinations, the best fit value resides in 
the phantom-like region, although a large area of quintessence-like region is also allowed. 

The IHDE model is mildly favored over the ACDM model according to the BE criterion, 
the evidence is slightly stronger than the HDE model case, but not yet sufficient for drawing 
strong conclusions. 

V. CONCLUSION 

In this paper we firstly gave a brief review of holographic dark energy model for both 
the non-interacting case and interacting case. We introduced a new interacting term Q = 
3ct-£/pde and the non-interacting case could be viewed as the special case with a = 0. We 
derived the equations for the evolution of Qde and H(z), and illustrated the dynamical 
behavior of these models by chosen some representative values of the parameters c and a. 
The condition for the model to have "Big Rip" is determined. 

Secondly, we utilize several data sets from the resent observations to constrain the models. 
Our data sets consist of the selected 182 high-quality type la supernovae, the baryon acoustic 
oscillation measurement from the Sloan Digital Sky Survey, the latest X-ray gas mass fraction 
data from Chandra observations, 27 GRB samples generated with E pea ^—E 1 correlation, and 
the CMB shift parameter from WMAP three years result. We used the MCMC technique 
to simulate the posterior probability of the model parameters. The best-fits for the three 
data sets are given in Table 1 for the HDE model, and Table 2 for the IHDE model. We 
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also give the probability distribution of the parameters and the contour maps for the HDE 
and IHDE models. 

Next, we utilize the Bayesian evidence (BE) as a model selection criterion to compare 
the holographic models with ACDM model for the three data sets. The BE is particularly 
appropriate for comparing models with different number of parameters. Both the HDE and 
the IHDE model are mildly favored by the current observational data set, although the 
evidence is weak. For both the HDE and IHDE models, the data favors "quintom" behavior 
slightly, i.e. the dark energy initially has w e ff > —1, but eventually crossing the phantom 
dividing line, and the model ends with a "Big Rip" . However, quintessence-like behavior is 
also still allowed with the present data. 

In brief, we conclude that according to the combined measurements data, the holographic 
dark energy model, especially the interacting holographic dark energy model is mildly favored 
by the observations, and for the best fit model the equation of state for both the HDE and 
IHDE crosses —1, for which the Universe ends up in a Big Rip. 
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